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Abstract 

We show that the complex phase diagram of high T c superconductors can be 
deduced from a simple symmetry principle, a SO (5) symmetry which unifies 
antiferromagnetism with d wave superconductivity. We derive the approxi- 
mate SO (5) symmetry from the microscopic Hamiltonian and show further- 
more that this symmetry becomes exact under the renormalization group flow 
towards a bicritical point. With the help of this symmetry, we construct a 
SO (5) quantum nonlinear a model to describe the effective low energy degrees 
of freedom of the high T c superconductors, and use it to deduce the phase dia- 
gram and the nature of the low lying collective excitations of the system. We 
argue that this model naturally explains the basic phenomenology of the high 
T c superconductors from the insulating to the underdoped and the optimally 
doped region. 
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In search of a low energy effective theory The high T c superconductors are among 
the most complex systems ever studied in condensed matter physics. In order to reach to 
the core of the physics responsible for this novel phenomenon, we obviously have to cross 
many different energy scales, and integrate away irrelevant degrees of freedom. Earlier in 
the high T c development, Anderson Zhang and Rice [Q] argued rather successfully that 
a good starting point for modeling the strong correlation effects in the oxides should be the 
Hubbard, or the t — J model close to half filling. Unfortunately, both of the these models 
are still too difficult treat. Their exact solutions in one dimension shed very little light on 
how to solve these models in higher dimensions. 

Even within the much simplified Hubbard or t — J model close to half filling, there are 
two different energy scales in the problem. There is a piece of "high energy physics" , which 
is responsible for forming local singlet pairs. This piece of physics is not very subtle, and 
basically originates from the J term in the Hamiltonian. We can associate a temperature, 
Tmf with this energy scale, which can be between zero and 1000-fT, depending on the filling 
factor. For conventional BCS superconductors, the energy scale of the pair formation is the 
same as the true phase transition temperature into a superconducting state. However, for 
the oxides, this piece of "high energy physics" does not tell us about the true ground state 
of the system. 

Back in 1987, Lee and Read asked a thought provoking question, "Why is the T c 
so low?". If the singlet formation is the same as superconducting transition, the natural 
transition temperature would be the same as Tmf-, much higher than the observed T c . The 
reason that the T c is low is because there is another "low energy physics", valid below T MF , 
which governs the fate of the singlet pair. The singlet pairs can form a spatially order 
state, leading to a antiferromagnetic ground state, they can condense, to form a d wave 
superconductor, they can arrange themselves into a spatially non-uniform state and lead 
to phase seperation Q], they can form a spatially homogenous mixed state of coexisting 
antiferromagnetic and superconducting order, a "spin bag phase" || , and finally, the effects 
of quantum fluctuation can be strong enough to destroy any form of order parameter, leading 
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to a "resonating valence bond (RVB) phase" jL|. In order to sort out the different competing 
ground states, we need a effective Hamiltonian which describes this "low energy physics" 
sector of the t — J model, valid below Tmf- It should be a simple Hamiltonian which can be 
solved analytically, and yet, the complexity of the possible phases should emerge from the 
unity of the model and the simplicity of its mathematical solution. 

What is the principle which can guide us in the search of such a universal low energy 
Hamiltonian? The different competing orders are not separated by any obvious energy scales 
in the remaining problem. Therefore, it is not obvious how one can systematically apply the 
renormalization group idea to integrate out the irrelevant degrees of freedom. Fortunately, 
whenever Nature poses us a extraordinary difficult challenge involving strongly interacting 
low energy degrees of freedom, it often times gives away its secret in the symmetry of the 
problem. Symmetry principle, exact or approximate, served us very well in our quest of the 
laws of nature, and it appears to be the only principle which can help us to solve the highly 
complex high T c problem. The main idea is to identify the symmetries of the microscopic 
Hamiltonian, such as the t — J model, and use the symmetry to constrain the possible form 
of the low energy effective Hamiltonian. 

What symmetries should we use, and what symmetries remain to be discovered? Well, we 
know the low energy effective Hamiltonian of a antiferromagnet, it is the 50(3) quantum 
nonlinear a model. Chakravaty, Halperin and Nelson |J showed convincingly that the 
physical properties of the undoped oxides are well described by this model. More recently, 
Pines and coworkers , Chubukov, Sachdev and Ye || argued that this model is applicable 
in the underdoped regime as well. We also know the low energy effective Hamiltonian of a 
superconductor, it is described by a £7(1) quantum nonlinear a model, or sometimes called 
the XY model. Doniach and Inui |TU| attempted to describe the metal insulator transition 
in terms of a quantum XY model, and recently, Emery and Kivelson |]TT[ presented evidence 



that the superconducting transition on the underdoped side of the oxides can be described 
by a renormalized classical XY model. Both the £0(3) spin rotation and the U(l) charge 
rotation symmetries are the obvious symmetries of the microscopic t — J model, and should 



certainly be used as constraining symmetries of the new, unified low energy theory 

The simplest way to construct a unified theory of antiferromagnetism and superconduc- 
tivity is to introduce a concept which I call a superspin^. It is a five dimensional vector, 
n a = (n>i, n2, ns, 714, ng). The first and the fifth components are the superconducting com- 
ponents of the superspin, identified with the two d wave superconducting order parameters 

0,1 



m = A f + A , n 5 = -z(A f - A) 

At = X)#(p) c It c -Pi ' g(p) = cosp x - cospy (1) 
p 



The remaining three components are the spin components of the superspin, identified with 
the antiferromagnetic order parameter 

n 2 = C p+Q,i a ij C Pj ' n 3 = X! C p+Q,i <T ^' C P.i ' n4 = C p+Q,i a ij C Pd (^) 
P P V 

where cx a, s are the Pauli spin matrices and Q = (n, 7T, 7r) is the antiferromagnetic ordering 
vector. The concept of the superspin is similar to that of the pseudospin |i3|JT4| , with a 



crucial difference that pseudospin is really "pseudo" in the sense that it has no real spin 
component. There is a 5*0(3) spin symmetry acting on the (ri2, ^3, ^4) subspace, with the 
total spin S a being the generator of the rotation, and there is a U(l) charge symmetry acting 
on the (711,71.5) subspace, with the total charge Q being the rotation generator, 



v 



S a = Y, c U<j c P,i > Q=\{N-M) (3) 



Here N is the number of electrons and M is the number of lattice sites. Of course, the concept 
of the superspin would only be useful if we can enlarge the known 50(3) x U(l) symmetry 
group, to include orthogonal transformations which can rotate the antiferromagnetic order 
parameters into the superconducting ones. The minimal group to accomplish this would be 



I would like to thank Prof. R.B. Laughlin for a stimulating discussion which helped me to form 
this concept. 
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a SO (5) symmetry group. If such orthogonal transformations exists, and if they commute 
with the microscopic Hamiltonian, then the concept of the superspin would be a very useful 
one. Since the orthogonal transformations preserve the magnitude of the vector, we can 
impose the constraint 

< = 1 (4) 

and really think of the superspin as a vector of a fixed magnitude in a five dimensional space, 
where the antiferromagnetic and superconducting order appear as three and two dimensional 
projections of this five dimensional vector. 

But we need 10 different "Euler angles" to describe a general rotation in a 5 dimensional 
space, mathematically speaking, the 50(5) group has 10 symmetry generators. (In general, 
the SO(n) group has n(n — l)/2 symmetry generators). The known 5*0(3) x U(l) symmetry 
supplies us with 4 such generators, where are the remaining 6 generators? 

Enter the 7r operators and the 50(5) algebra Fortunately, the remaining 6 symme- 
try generator have already been discovered by Demler and Zhang flT5| |, in their recent work 



on the theory of the resonant neutron scattering of the high T c superconductors. These 
operators carry charge 2, spin 1, and a total momentum (7r, tt, it). Here we will call them ir 
operators, and they are defined as follows 

*l = Y.9(pH +Q ^ a ° y h c -P,j > = (5) 

p 

and we see that there are obviously six of them. 

Let us recall that the 50 (5) Lie algebra has ten generators L ab = —L ba satisfying the 
following commutation relations 

[L a bi L c d] = i($acLbd ~ ^ad^bc ~ hcL a d + ^bdL a c) (6) 

Let us identify the ten operators (Q, S a , ir a , njj with the generators of the 5*0(5), L a b, in 
the following way, 
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J ab 







7r| + Tl x 



\ 



(7) 





-s z 

Q -z(tt| - vr^) -z(ttJ - 7r„) -z(tt] - tt z ) y 

where the value of the matrix elements on the upper right triangle are determined by anti- 
symmetry, L ab = —L ba . Using the fermionic representations of the (Q, S a , n a , tt^) in equa- 
tions (§), (H), and the above identification, we can check explicitly that the 5*0(5) commuta- 
tion relations (0) are indeed satisfied provided that we make the continuum approximation 



(8) 



This approximation is valid when we consider the long wave length limit, where the only 
important information about the internal structure of a composite operator is its symmetry, 
but its detailed shape. Since g 2 (p) has a s wave symmetry, we approximate it by its average 
value. 

With the aid of the 5*0(5) algebra, the concept of the superspin introduced earlier can 
be defined precisely. The superspin n a is defined as a vector representation of the 50(5) Lie 
algebra, satisfying the following commutation relation 



[L a bi n c 



-i8 bc n a + i8 ac n b 



(9) 



If we substitute for L ab and n c the microscopic definitions in terms of the electron operators 
(H), (||) and (0), we find that they are indeed satisfied in the continuum approximation (||). 
In particular, we see that the 7r operators rotate the antiferromagnetic order parameter into 
superconducting order parameter and vice versa. Therefore, the six 7r operators |I5J recently 
introduced by Demler and Zhang are exactly the six missing symmetry generators we are 
searching for. 

With the basic low energy variables assembled and their kinematic relations fully speci- 
fied, we are now in a position to address the dynamic properties, and discuss the commuta- 



tion relations between the symmetry generators L a b and the microscopic t — J Hamiltonian. 
The total spin operators L 2 3, L 2 4, L 34 and the total charge operator L 15 certainly commute 
exactly with the t — J Hamiltonian. What about the remaining six generator? 

Demler and Zhang |15] have shown that although they do not commute exactly with the 
Hamiltonian, they are approximate eigenoperators of the t — J Hamiltonian, in the sense 
that 

[H,ni}=u 7ii (10) 

where ujq = J(l — n)/2 — 2fi. The way to obtain this relation is by the standard t matrix 
approximation, similar to the random-phase-approximation for particle-hole operators. The 
reason that this relation hold is rather simple. The irj^ operator is a particle particle opera- 
tor, describing a pair of electrons with a center of mass momentum Q = (jr, n, it). However, 
for tight binding models with nearest neighbor hopping, the two particle continuum vanishes 
at Q. Therefore, a repulsive interaction, such as J in the triplet channel, naturally leads 
to a antibound state, uq in the above equation is nothing but the energy eigenvalue of this 
antibound state. More sophisticated calculations beyond the t matrix approximation have 
also been carried out recently P^JT7[ ]. Equation ( |TUD is a fundamental equation underlying 
the present work, it provides a bridge from the microscopic t — J Hamiltonian to the effective 
nonlinear a model theory. One can also define a tc' operator, using g(p) = cosp x + cosp y in 
equation (|5|). Such a operator would transform antiferromagnetism to s wave superconduc- 



tivity. However, it was found |T5| that they are in general not eigenoperators. Therfore, our 
theory naturally predicts d pairing symmetry. 

There are two reasons why the SO (5) symmetry is only a approximate one. One is the 
approximation used to derived equation flTop. Therefore, this relation should be tested in 
exact diagonalization studies of the Hubbard or the t — J model. This can be done by 
calculating the dynamical correlation functions of the 7r operators and identifying isolated 
poles in their spectral function. Fortunately, such a study has been carried out recently by 
Meixner and Hanke |18|| , there is indeed a sharp collective mode in the ir channel. The second 



reason why the SO (5) symmetry is broken explicitly is that even if ([H]) holds exactly, the 
Hamiltonian does not commute with all generators of 50(5). However, this source of explicit 
symmetry breaking is easy to handle, and has interesting physical consequences which we 
will explore throughout this paper. The eigenoperator relation (|H]) is almost as good as the 
vanishing of the commutator. The reason is that the Casimir operator C = J2 a <b L 2 a b °f the 
SO (5) algebra commutes with the Hamiltonian, and therefore, we can still classify all the 
eigenstates of the Hamiltonian according to their SO (5) quantum numbers. This also means 
that the superspin vector multiplet n a is not mixed into other higher tensorial representations 
of the 50(5) under time evolution, a crucial fact which enables us to construct a low energy 
effective Hamiltonian containing the superspin degrees of freedom only. The situation is 
very analogous to the problem of a spin precessing in a magnetic field. While the full spin 
rotation symmetry is broken by the magnetic field, the magnitude of the spin is preserved. 

So far I have tried to present the line of arguments according to the logical reasoning in 
the high T c problem. However, at this stage, it would be appropriate to reveal the historic 
origin and the source of inspiration which lead to the present work. The present work is a 
generalization of the concept of the exact 50(4) symmetry |nj] of the Hubbard model and its 
application to the negative U Hubbard model []HpO| ,ZT|. Using the rj operators constructed 



earlier by Yang [ 22| , Yang and Zhang [ T9j pointed out that the Hubbard model has in 



additional to the usual 50(3) spin rotation symmetry also a 50(3) pseudospin symmetry 
generated by the following operator, 

= c I+QT c -pI ' r l = (v j V , Vo = Q (11) 
p 

While this symmetry is valid both for the positive and the negative U Hubbard model, the 
simplest vector multiplet contains s wave superconducting and charge density wave order 
parameter 0,|2O[, 

A i = E4Ai > ^o = J2 c l +Q ,Aa , A. = (Aj)t (12) 

p pa 

and is therefore only useful for the negative U model. The simplest generalization of the 



pseudospin symmetry to the positive U or the t — J model is the concept of the super- 
spin introduced in this work, the precise correspondence between these two symmetries are 
summarized in table I. 

Construction of the 5*0(5) quantum nonlinear a model With the basic kinematic 
and dynamic constraints assembled, we are in a good position to construct a low energy 
effective Hamiltonian. But before we do so, a critical reader may be concerned with the fact 
that the SO (5) symmetry is only approximately valid. How do we address this question? 
Well, it may be useful to look at a historical example where a approximate symmetry was 
used solve a important physics problem. The quantum chromodynamics (QCD) is a model of 
strong interaction. Unlike the quantum electrodynamics, it can not be solved perturbatively. 
However, in the limit of vanishing quark masses, the QCD Lagrangian has a SU(2) x SU(2) 
symmetry, with conserved vector and axial currents. In real world, the quarks have small 
(compared to the natural QCD energy scale) but finite masses, so that the axial current 
is only approximately conserved. The idea behind the hypothesis of partially conserved 
axial current (PCAC) is to construct a effective nonlinear a model assuming that the axial 
currents were conserved, and then add small explicit symmetry breaking terms to reflect 
the effects of finite quark masses. The nonperturbative, low energy dynamics of the QCD 
is essentially understood in this way. Therefore, we will follow the same strategy here to 
extract the low energy content of the t — J model. Along the way, we shall encounter a 
pleasant surprise, although the SO (5) symmetry that we start from is approximate, near a 
critical point, it becomes exact under the renormalization group flow. 

Let us first assume that the SO (5) symmetry were a exact symmetry, and ask what is 
the form of the low energy Hamiltonian. As we argued before, below T M f, the system has 
a tendency to form local singlets, but it has not fully decided whether to become antiferro- 
magnetic or superconducting in its ground state. Therefore, we can think of T MF as a mean 
field transition below which the superspin acquires a fixed magnitude, leaving its orientation 
as a low energy degrees of freedom. Close to Tmf, the anisotropy terms in the superspin 
space are not very important, and this transition can be simply described by a standard 
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Landau free energy functional 

F = a\n\ 2 + b\n\ 4 (13) 

T MF is the temperature below which the coefficient of the quadratic term a = a'(T — T MF ) 
changes sign. 

Below T MF , the magnitude of the superspin is fixed, we can always rescale it to satisfy 
the constraint (|J). The ten generators of 5*0(5) describe the orthogonal rotation of the su- 
perspin, leaving its magnitude fixed. Having either antiferromagnetism or superconductivity 
corresponds to a fixed direction of the superspin, thus breaking the 5*0(5) symmetry sponta- 
neously. The low energy dynamics of such a system is determined completely in terms of the 
Goldstone bosons and their nonlinear interactions as specified by the 50(5) symmetry. The 
kinetic energy of the system is simply that of a 50(5) rigid rotor, given by ^ J2 a <b L 2 ab {x). 
The quantity x is nothing but the moment of inertia of the rigid rotor. If we assume the 
exact 50(5) invariance, and rescale the magnitude of the superspin to satisfy the constraint 
(|j), there are no polynomial terms in the potential energy, and in the long wave length limit, 
we can perform a gradient expansion to obtain a term ^ J2 a {dkn a {x)) 2 to the leading order. 
Therefore, the resulting Hamiltonian density is given by 

H s = E L 2 ab (x) + P - £(<W*)) 2 (14) 

Z X a<b Z a 

H s is a unique low energy Hamiltonian constrained by the 50(5) symmetry. However, 
this symmetry is not exact, therefore, we should allow some weak symmetry breaking pertur- 
bations. Of course one should break the 50(5) symmetry in such a way that the subgroup 
50(3) x U(l) of the spin and charge symmetry is still exact. Physically, this corresponds 
to going from a symmetric rotor of only one moment of inertia to a asymmetric rotor with 
different moments of inertia. The asymmetric Hamiltonian is given by 

H a = J2 ^-L ab (xf + £ ^(v k ab (x)) 2 + V(n) (15) 

a<b ^Xab a<b * 

where xis = Xc X23 = X24 = X34 = Xs and Xi(2,3,4) = X(2,3,4)5 = Xn are the charge, spin and 
the newly introduced 'V susceptibility, and we have also introduced a generalized velocity 
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field v* b = n a dkn,b — ribdkn a and the corresponding stiffness in the charge (pis = p c ), spin 
(P23 = P24 = P34 = Ps) and the ix (pi( 2 ,3,4) = P(2,3,4)5 = Ptt) channel. In the presence of 
explicit symmetry breaking, we also allow a quadratic symmetry breaking term of the form 

V(n) = - 9 -(nj + nj + nl) (16) 

H s and H a describe the system at half filling only. In order to go away from half 
filling, one simply has to add a chemical potential term —2pQ = — 2pL^, to obtain 7i s = 
H s — 2pLi5 and 7i a = H a — 1\xL\§ respectively. This Hamiltonian is to be quantized using 
the commutation relations (||) and @, and a complete set of equations of motion can be 
determined. 

Given the Hamiltonian 7i a , one can simply perform a Legendre transformation and obtain 
the corresponding Lagrangian. After a simple Wick rotation t — > ir, the Lagrangian density 
becomes 

£ a = E ^ ab (x) 2 + E f»)) 2 + V(n) (17) 

a<b a<b 

where 

u a b = n a (d T n b - iB bc n c ) - (a — >• b) (18) 

is the angular velocity. Here we have introduced a set of generalized external potentials B ab 
coupling to L ab . In our current problem, the only nonvanishing component is B i5 = —B 5 i = 
2/i. The partition function of the system is given by 

Z = J [dn a }5(n 2 - l) e ~ lo dT I (19) 

It is now appropriate to define the constants appearing in the model. Xs and \c are 
the familiar uniform spin susceptibility and charge compressibility of the system, while p s 
and p c are the spin stiffness and charge stiffness respectively (p s should not be confused 
with the superfluid density which sometimes also uses this notation). Xn and p^ are the two 
constants newly introduced in this work. They describe the time and length scale over which 
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a antiferromagnetic region can be converted to a superconducting region and vice versa. In 
order for the SO (5) symmetry to be approximately valid, one would require these constants 
are close in value. The last remaining constant of our model is the anisotropy constant g, 
which selects either a "easy plane" in the superconducting space (711,715), or a "easy sphere" 
in the antiferromagnetic space (n 2 ,n 3 ,n 4 ), depending on the sign of g. As argued before, 
/i = defines the model at half filling, where we know the system is antiferromagnetic. 
Therefore, we fix g > to match this fact. At this moment, the values of the constants 
(especially the sign of g) of the model are not deduced from any first principle calculations. 
However, once these values are fixed at half filling, we will not allow them to vary in a 
arbitrary fashion. As we shall see in the next few sections, the richness of phase diagram 
comes entirely from the variation of the chemical potential. 

Origin of superconductivity At half filling, we choose g > so that the superspin 
prefers to lie in a "easy sphere" of (712,713,714). Away from half filling, the only new term 
appearing in the Hamiltonian is just the chemical potential term in 7i a . When we investigate 
the Lagrangian C a , we observe that the chemical potential term appear as a "gauge coupling" 
in the imaginary time direction. A constant chemical potential term is a "pure gauge" and 
could therefore be "gauged away". One's naive expectation would be that such a term has 
no dynamic consequences. However, one has to pay a price for doing so. Although a constant 
chemical potential term can be gauged away in the bulk, it reappears as a twisted boundary 
condition in the imaginary time direction, and could have non-trivial consequences. 

It is more direct to investigate the Lagrangian (pj]) with the periodic boundary condition 
in the imaginary time direction. Because of the periodic boundary condition, the classical 
path extremizing the path integral are the static solutions. For the static solutions, we see 
that the only nonvanishing contribution of the kinetic energy is the chemical potential term, 
which gives a effective potential energy 

V eff = V(n) - M^(n? + nt)[ Xc (nl + nj) + Xir (n* + nj + n\)} (20) 

Let us first take the 50(5) symmetric situation where x-k = Xc — X? i n this case the terms 
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in the square bracket reduces to x- If we were dealing with a abelian XY model, such a 
term still does not have any dynamical consequences. Since n\ + n\ = 1 in such a case, the 
above term reduces to a number and gives only a trivial shift of the ground state energy, 
and can not lead to any nontrivial phase transitions. However, if the XY symmetry is 
embedded into a higher symmetry group, as is done in the present case, this term has a 
profound dynamic consequence. While the g term with g > favors a antiferromagnetic 
"easy sphere" (n 2 , n^, 714), the chemical potential term favors a superconducting "easy plane" 
(ni, ns). This competition leads to a first order (for more detailed discussion of the order of 
the transition, see a later section) phase transition when 

1 



v = Vc = -yg/x (21) 

For n < /i c , the ground state is antiferromagnetic, and the superconducting state has a 
finite energy. This energy decreases gradually as one increases /1, until a level crossing at 
fi c . For /i > /i c , the situation is reversed, the ground state is superconducting, and the 
antiferromagnetic state has a finite energy. 

For Xn 7^ Xc, the situation is a bit more complex. In the parameter regime Xc > Xtt, 
there is always a first order direct transition from antiferromagnetism to superconductivity 



at fi = /i c = \\j g/Xc- On the other hand, in the parameter regime Xc < X-n- < 2Xc, and for 

— < (2fi) 2 < § _ 1 (22) 

XlT ^ Xc Xtt/ ^ 

there exists a intermediate "spin bag" phase with coexisting antiferromagnetic and su- 
perconducting order. Schrieffer, Wen and Zhang || describe such a phase in terms of pairing 
the eigenstates of the antiferromagnetic background. It is amusing to observe that when 
expressed back in terms of the orginal electron operators, their order parameter is a mixture 
of the antiferromagnetic, d wave superconducting order parameter and the 7r operators. 

The above discussion based on the Lagrangian may appear rather formal for some read- 
ers, a intuitive physical example would be useful. Consider a antiferromagnet with a easy 
axis anisotropy (say the z axis). Below the Neel transition, the Neel vector prefers to point 
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along the easy axis. A uniform magnetic field pointed along the easy axis creates a easy xy 
plane. At a critical value of the magnetic field, there is a "spin flop" transition where the 
Neel vector changes its orientation from the z axis to the xy plane. Once the underlying 
50(5) symmetry is revealed, the physics of the high T c superconductivity is as simple as 
the "spin flop" transition. We see that equation flTgp describes a precession of the superspin 
about a "fictitious magnetic field", namely the chemical potential. The easy axis of the 
antiferromagnet translates into the "easy sphere" {11,2,113,114) while the "easy plane" of the 
antiferromagnet translates into the "easy superconducting plane" {ni,n§). The transition 
from a antiferromagnetic ground state to a superconducting one is that of a "superspin flop 
transition" . In fact the phase diagram of a easy axis antiferromagnet in the temperature 
and magnetic field plane is very similar to that of a high T c superconductor in the temper- 
ature and chemical potential plane. The complete analogy between these two systems are 
summarized in the following table. 





t — J model 


negative U Hubbard 
model 


ant if err omagnet 
in a B field 


symmetry 


SO{5) 


50 (4) 


S0{3) 


symmetry generators 




S a ,Q,ri,rjt 


S a 


order parameter 


superspin 


pseudospin 


Neel vector 


symmetry breaking 


n 


n 


B 


phase transition 


superspin flop 
(from SDW to d wave 
superconductivity) 


pseudospin flop 
(from CDW to s wave 
superconductivity) 


spin flop 
(from easy axis to 
easy plane) 


collective modes 


1 massless phase and 
3 massive SDW modes 


1 massless phase and 
1 massive CDW mode 


1 massless XY and 
1 massive Z mode 



This analogy also helps us to understand the origin of the superspin flop transition in 
the Hamiltonian formulation of the problem. In the Hamiltonian formulation, the chemical 
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potential term appear as a coupling to the symmetry generator L15, the number operator Q, 
but not to the superspin n a directly It is not immediately obvious why this term would select 
any particular direction in the superspin space. Of course, the same problem occurs in the 
antiferromagnetic analogy, where the uniform magnetic field couples to the total spin vector, 
not the Neel vector. If one tries to visualize the the Neel Hamiltonian in a semi-classical 
fashion, one has to remember that there is a important constraint that the total spin vector 
is orthogonal to the Neel vector. The externally applied uniform magnetic field leads to a 
finite total spin in the z direction. In order to satisfy the orthogonality condition, the Neel 
vector therefore has to lie in the plane orthogonal to z. In the present model, (restricting to 
the SO (5) symmetric case to simplify the algebra), the orthogonality constraint takes the 
form 

e abcde n c L de = (23) 

which can be simply proved by expressing the angular momenta in terms of the angular 
velocities L a b = x( n a^& — Wf,n a ). The chemical potential term leads to doping, or prefer a 
finite value L 15 . We immediately see that the constraint equation gives n 2 = n 3 = = if 
other L a b generator have no ground state expectation value, which is the case since x > 0- 
The question of the origin of superconductivity is not one question, but rather two related 
questions separated by a energy scale. The "high energy" mechanism leads to binding of 
electrons into singlet pairs. The origin of this binding is rather obvious in the t — J model, 
since the J term favors electrons on near neighbor sites to be antiparallel. The problem 
is that this same J could also lead to antiferromangetism. In our new language, this pair 
binding gives rise to a finite magnitude of the superspin, leaving its orientation still unfixed. 
This is in marked contrast to the celebrated Bardeen-Cooper-Schrieffer theory of conven- 
tional superconductivity where pair binding is equated with the onset of superconductivity, 
and also different from the large phase fluctuation picture [JD]] where the only uncertainty 
of the order parameter is its phase. In addition to the high energy pair binding mechanism, 
a "low energy" mechanism selects a orientation of the superspin and distinguishes antiferro- 
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magnetism from superconductivity. The selection of the different possible ground states are 
described by 50(5) nonlinear a model and the mechanism for favoring superconductivity is 
the "superspin flop" mechanism discussed above. We see explicitly that superconductivity 
is a inevitable consequence of the 50(5) symmetry and antiferromagnetism at half filling. 
From this point of view, a antiferromagnet state can in some sense be thought of as a solid 
formed by Cooper pairs, and the "superspin flop" transition is a first order melting transition 
from the solid into a superfluid of the Cooper pairs. In this frame work, the "spin fluctua- 
tion exchange" calculation should be interpreted as a mean field theory of Tmf- Such 
calculations are important parts of our understanding, and lead to the glorious prediction of 
d superconductivity, but they are not complete since they do not address the actual phase 
transition and the competition between antiferromagnetism and superconductivity. The two 
different mechanisms compliment each other in their respective energy regime and together 
form a complete picture of the origin of superconductivity in the oxides. 

Theory of collective modes and their nonlinear interactions The superspin model 
gives a natural description of the collective modes in the ordered phase. Close to the tran- 
sition between the antiferromagnetic and the superconducting ground states, the ordered 
phases are not conventional, but their low energy excitations reflect the competition between 
the two different kinds of ordering. These low energy excitations can be studied systemat- 
ically by applying symmetry principles and they have profound experimental consequences 
which we shall discuss in this section. 

Before performing concrete calculations lets us first address the qualitative features of 
the low energy excitations using symmetry principles. Let us first suppose that the 50(5) 
symmetry were exact. In this case, if the superspin has a fixed orientation in its ground 
state, this state breaks the 50(5) symmetry spontaneously. Spontaneous breaking of a 
continuous symmetry naturally leads to Goldstone bosons. The number of the Goldstone 
bosons is the number of the broken symmetry generators. A fixed direction of the superspin 
leaves the 50(4) subgroup of the 50(5) unbroken. The 50(5) group has 10 generators, 
while the 50(4) group has 6 generators. Therefore, we would expect to have 4 Goldstone 
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bosons, corresponding to the 4 broken symmetry generators. 

What is the nature of these 4 Goldstone bosons? Suppose that the superspin vector lies 
in the superconducting plane, say along the n\ direction. The mode corresponding to the 
broken rotation generator Li 5 is just the usual Goldstone mode, describing the phase degrees 
of freedom of the superconducting order parameter. The three other modes corresponds to 
the broken generators Li2,Li3 and L14. These modes form a triplet representation of the 
unbroken 5*0(3) spin symmetry. On the other hand, if the superspin vector lies in the 
antiferromagnetic sphere, say along the 112 direction, the usual spin waves corresponds to 
the broken generators L23 and L24. They do not form a representation of the spin 5*0(3) 
group, since this symmetry is spontaneously broken. In addition to the usual spin waves, 
the current theory predicts two additional Goldstone modes, corresponding to the broken 
generators L 2 i and L25. They form a doublet representation of the charge U(l) symmetry 
group. 

This symmetry analysis make it clear what happens if the 50(5) symmetry is explicitly 
broken to 50(3) x U(l). In the presence of a quadratic symmetry breaking term g e ff = 
g — x(2/i) 2 , the spin triplet Goldstone mode of the superconducting state would become 
massive if g e ff < and the charge doublet mode of the antiferromagnetic state would 
become massive if g e ff > 0. The usual spin wave modes and the phase mode remain 
massless because of the 50(3) x £7(1) symmetry. Goldstone bosons which become massive 
because of explicit symmetry breaking terms are sometimes called pseudo Goldstone bosons. 

This general symmetry analysis can be easily checked by a explicit calculation starting 
from the 7i a . The equations of motion derived from this Hamiltonian takes the form 

Lab{x) = ~ J2(Xac ~ Xbc){ L ac(x), L bc (x)} + ^(p aC ~ Pb^V^{x)^{x) 
Z c c 

+dk ^2(Pa C n b n c v^ c (x) - p bc n a n c vl c {x)) + ^2(B ac L bc (x) - B bc L ac (x)) 

c c 

+9 ( n a$b a n a - n b 5 aa n a ) (24) 

a=2,3,4 

Lab(x) = Xab^ab (25) 

where u ab is given by the Wick rotation of equation fli"8D . The first term on the right hand 
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side describe the local motion of a five dimensional asymmetrical top, a simple generalization 
from the equations of motion of a asymmetrical top in classical rigid body mechanics. The 
second term arises from the asymmetrical coupling of the neighboring rotors. The third term 
is the most relevant one for our discussion. It arises mathematically from the Schwinger term, 
proportional to the derivative of a delta function, in the local commutator between L a &(x) 
and v^ d (y), and expresses the anomaly content of the SO (5) current algebra. In the absence 
of the other terms on the right hand side, there are ten conserved charges, and this term 
describes the conserved currents associated with these charges. The fourth term simply 
describes the precessional motion of L a b(x) in the presence of a external field, and the last 
term describe the local torque due to the anisotropy potential. It is easy to see that for Q(x) 
and S a (x), only the third term is nonvanishing, expressing the continuity of these exactly 
conserved charges. 

This set of equations are the nonabelian generalization of the familiar Josephson equa- 
tions [^] in conventional superconductors. The first equation is the generalization of the 
Josephson current relation. If we take L a b(x) to be L 15 (x) in the first equation, we find that 
the only surviving term is the third one, which reduces to the standard Josephson current 
expression J — p c W in the isotropic limit. The second equation is the generalization of the 
Josephson acceleration equation. If we take L a b(x) to be L\^{x) in this equation, we recover 
the familiar ac Josephson equation fi = d6/dt. This set of equations are also generalization 
of the Landau-Lifshitz-Neel equation of antiferromagnets p4 |. This statement can be easily 
checked if we take L ab (x) to be L 2 3(x), L 2 ^{x) and L^(x) in the above two equations. 

This set of equations are the fundamental equations of the high T c superconductors. 
They unify the treatment of antiferromagnetim with superconductivity and describe the lin- 
ear spectrum and the nonlinear interactions of the collective modes in the antiferromagnetic 
and the superconducting phases. The new physics contained in these equation can be easily 
analyzed in the linearized approximation. For g e ff > 0, the ground state is antiferromag- 
netic, we can linearize the equations around the Neel vector, say n 2 , and we obtain 
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Xsd 2 t n a = p s d 2 k n a , a = 3, 4 ; x^ 2 t n x = p^d 2 k n x - (g - (2p) 2 Xn)n x , A = 1, 5 (26) 

The first equation describe the two spin wave modes of the antiferromagnet, while the 
second equation predicts a new, massive doublet pairing mode. This mode is the precursor 
of superconductivity in the antiferromagnetic phase. This mode also makes it plausible that 
one can think of a antiferromagnet as a quantum solid of Cooper pairs. A quantum solid 
has two types of excitations, the gapless phonon modes, and a gapful mode corresponding to 
extracting a atom from its position |23| . The pairing mode described in the above equation 
roughly corresponds to the second type of excitation of a solid, if we identify the Cooper 
pair as a "atom" of the solid. We see that as the chemical potential is increased, this pairing 
mode lowers its energy, the quantum solid becomes softer, and eventually the solid melts 
completely to form a superfluid of Cooper pairs. 

For g e jf < 0, the ground state is superconducting, and we can linearize the equations 
around the n\ vector to obtain 

Xcd 2 n 5 = p c d 2 k n 5 ; x^d 2 n a = p^d 2 k n a - ((2 / u)\ 7r - g)n a , a = 2, 3, 4 (27) 

The first equation describe the usual Goldstone mode (sound mode) of a superconductor, 
while the second equation describe the triplet of massive magnetic modes predicted by Dern- 



ier and Zhang ||15|| . This mode is a precursor of antiferromagnetism in the superconducting 
phase. It can be roughly thought of as the roton mode in a superfluid, since they both reflect 
the "diagonal short range order" in a superfluid. As the chemical potential is decreased, this 
mode lowers its energy and eventually the superfluid "solidifies" to form a antiferromagnet. 

Both the doublet pairing mode in the antiferromagnetic phase and the triplet magnetic 
mode in the superconducting phase owe their existance to the kinetic energy terms of the ir 
operators. Therefore, we shall call them the tt doublet and the n triplet modes respectively. 

Both classes of new collective modes have important experimental consequences. The 
7i triplet mode has been used by Demler and Zhang []15] to explain the recent resonant 
neutron scattering experiments on Y BCO superconductors. It was found that below the 
superconducting transition temperature, a collective excitation peak appears in the neutron 
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scattering cross section. The energies of the mode range from 25meV to AlmeV [p5|-|2S 
depending on doping level, and the scatterings are observed in the spin triplet channel, at 
the commensurate momentum (tc, it, it), in exact agreement with the quantum numbers found 



theoretically. These modes vanish exactly at T c p7| , p8| , consistent with the interpretation 
that they are the pseudo-Goldstone modes associated with the U(l) symmetry breaking. 
The precise correlation between the energy of the neutron peak and the superconducting 
transition temperature will be discussed in a later section. 

The 7r doublet mode of a doped antiferromagnetic state is much harder to observe ex- 
perimentally. But they have indirect consequences which can be measured. In a puzzling 
experiment, Suzuki it et al. |30| observed a unusually long Josephson coupling length of up 



to 480A in the Y Ba2CujPij ' PrBa2Cu^Oi/YBa2Cu^O-j junction. This experiment could 
be interpreted as a signature of the small gap towards the tt doublet pairing excitation in 
the insulating material PrBa2Cu^OT. 

All above discussion assume that quantum fluctuations are not strong enough to destroy 
the long range order. In the opposite limit, when x^ 1 >> p, one should diagonalize the 
kinetic terms in 7i s and 7i a first and treat the potenial energy terms as perturbation. The 
kinetic terms can be diagonalized easily by using the representation theory of SO(n). They 
classified by a integer /, and the Casimir operator has the value of 1(1 + n — 2) in this 
representation. The ground state is a 5*0(5) singlet with I = 0, and the lowest energy 
excitation is a massive vector multiplet with 1 = 1. This multiplet is split by the anisotropy 
in the kinetic coefficients. The massive pairing doublet has a energy of ^ — h tt- while the 
massive magnetic triplet has a energy of — V 

The global phase diagram of high T c superconductors It is straightforward to 
determine the phase diagram of the 50(5) nonlinear a model and compare it with the 
experimentally observed phase diagram of high T c superconductors. The topology of the 
phase diagram takes different form depending on the strength of the quantum fluctuation. 
In the following, we will always be discussing the case of d = 3, unless otherwise stated. Our 
model does not have any disorder, therefore any disorder induced phase such as spin glass 

20 



is precluded in the current discussion. 

Let us first discuss the so called renormalized classical regime |J, where quantum fluc- 
tuations merely renormalize the coupling constants of the model, but they are not strong 
enough to destroy the various types of ordering. In this case, at half filling, fx = and g > 0, 
the superspin lies in the "easy sphere" (n 2 , n 3 , ra 4 ) and the ground state is antiferromagnetic. 
There is a phase transition at Tn to a paramagnetic phase as the temperature is raised, and 
this phase transition is in the universality class of the classical 50(3) model in three di- 
mensions. When the chemical potential is increased from its half filling value of fx = 0, the 
gap towards the 7r doublet mode decreases. Since it is easier for the superspin to fluctuate 
into other directions, the effective spin stiffness decreases and so does the Neel temperature. 
At the critical value fx = /i c , the superspin has no preferential "easy directions", and the 
system is described by a isotropic SO (5) nonlinear a model. This model has its own critical 
temperature Tb c , which is in general smaller than the Neel temperature, since the critical 
temperature of the SO(N) model scales inversely with N. If the fx — fi c line is reached when 
T < Tb c , the system is ordered and the superspin defined in this limiting procedure still 
lies within antiferromagnetic "easy sphere" . When fx is increased beyond the critical value, 
the superspin "flops" into the superconducting plane. Since the direction of the superspin 
changes abruptly at the transition, and since at the transition point, the correlation length 
in the sense of Josephson is finite, the "superspin flop" transition in general first order for 
T < T bc . Beyond the superspin flop transition, the system crosses from a 5*0(5) critical 
behavior over to a XY critical behavior, and one expects the superconducting transition 
temperature to increase with increasing chemical potential, since the mass of the 7r triplet 
mode increases. The topology of the phase diagram in this case is depicted in Fig. (la). 

The line of the first order superspin flop transition is in general not exactly vertical, since 
the susceptibilities and the stiffness in different superspin directions may not be the same. 
At zero temperature and for fx < fx c , Q — 0, the density of holes vanish. For fx > fi c , the 
density of holes is finite. The discontinuous jump of the density is given by 2x/i c . The reason 
for the hole density to vanish when fx < fx c is because we employed a superspin model with 
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a fixed magnitude. When one relaxes the model to a "soft superspin" one, the hole density 
can be finite in this region. However, the discontinuous jump in the density will remain in 
the "soft superspin" model. Because there is a discontinuous jump in the density, for hole 
concentrations in the "forbidden region" < Q < 2xfi c , the system phase separates into a 
hole rich and a hole deficient region, a possibility first pointed out by Emery, Kivelson and 
Lin @] in the context of the t — J model. The physics in the "forbidden density region" could 
be rather complicated, and interactions not included in our current model could give rise to 
various types of ordering, such as incommensurate order and "stripe phases" . However, we 
believe that this types of order are not directly related to the origin of superconductivity. 

The first order superspin flop line terminates at a critical point, from which two second 
order lines, one corresponding to the Neel transition and the other corresponding to the 
XY transition, emerge. For this reason, this critical point is called a bi-critical point. 
There are two relevant operators in the close vicinity of the bi-critical point, the reduced 
temperature t = (T — T^/T^ and g e ff- The scaling behavior close to this point is well 
studied in the literature |TI|-[33J . One of the most important result is that even if the model 
defined at the short length scales is not exactly SO(N) invariant, i.e. the susceptibility and 
stiffness in different directions may not be the same, the difference among them scales to 
zero under the renormalization group flow when the bi-critical point is reached. A simple 



version of this result is shown by Pelcovit and Nelson [32], who considers a SO(N) model 
where the the N — th component has a different stiffness than the others. They have shown 
that the difference operator has negative scaling dimension in the 2 + e expansion of the 
nonlinear a model, and is therefore irrelevant close to the critical point. A more general 
and very elegant result was derived by Friedan [|35|| , who considers a nonlinear a model in 
a arbitrary Riemanian manifold with metric g a b. Friedan showed that the fixed point of 
the renormalization group flow is a Einstein manifold with g ab = 27 r(d-2) -^ afe ' w here R a b is 
the Ricci tensor. In particular, if one starts with a distorted sphere S n , only the average 
curvature is a relevant variable, whereas the distortions preserving the average curvature are 
irrelevant. Therefore, if one has a microscopic model without the full spherical symmetry, 
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a perfect sphere with constant curvature is produced by the large scale fluctuations close 
to the critical point. This result provides a strong theoretical justification to describe the 
high T c superconductors by a SO (5) nonlinear a model: Even if the symmetry is only 
approximately valid in the microscopic models such as the Hubbard or the t — J model, close 
to the bi-critical point, where the most interesting transition from a antiferromagnetic state 
to a superconducting state occurs, this symmetry becomes exact. 

Close to the bi critical point, thermodynamic quantities obey exact scaling relations. 
As mentioned before, there are two relevant operators in the problem. Finite t gives a 
temperature correlation length according to the scaling law £t ~ t~ v . The other relevant 
parameter is the biquadratic symmetry breaking term g e ff, which gives a correlation length 
£g ~ 9eff Close to the bicritical point, the singular part of the free energy scales according 

f(t,g eff ) ~ t 2 - a f(Ag eff /t<t>) (28) 

where a is the free energy exponent of the 5*0(5) model and is a crossover exponent given 
by = v/vg. In equation (25), A is a non-universal constant and / is a universal scaling 



function of a single argument. Other singular thermodynamic quantities are be determined 
in a similar way. 

The crossover exponent <fi determines the way in which the two second order lines merge 
into the first order line. The universal scaling function / diverges at two arguments, x + > 
and x_ < 0. The two second order lines are determined by 

geff = {x + /A)1+ , g eff = (x-/A)t* (29) 

Explicit calculations of the exponents can be carried out within the large N approximation 
p6| . To the leading order, v — 1 and v g = 1/2, therefore <fi = 2. From this we conclude that 
the two second order lines merge into the first order line tangentially. (See Figure (la)). 

As we have seen before, there exist parameter regimes where antiferromagnetism and 
superconductivity can both coexist. In this case, the topology of the phase diagram has to 
be modified to include the immediate mixed phase, the "spin bag" phase ||. The "spin bag" 
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phase is conceptually similar to the conjectured supersolid phase in A H e . In this case, part or 
the entire first order superspin flop lines separates into two second order lines enclosing the 
"spin bag" phase. (See Figure (lb)). The four second order lines intersect at a tetracritical 
point T tc . 

So far all phase transitions discussed are classical phase transitions, in the sense that 
quantum mechanical fluctuations only renormalize the coupling constants in the problem. 
When quantum fluctuation is gradually increased, ordering in some phases can be destroyed. 
Since the SO (5) isotropic point has the largest quantum fluctuation, one would expect T^ c 
to be driven to zero first. In this case, the bicritical point is a quantum critical point where 
the antiferromagnetic state goes into the superconducting state via a direct second order 
phase transition. (See Figure (lc)). 

When the quantum fluctuations are increased further, a entire region between the antifer- 
romagnetic phase and the superconducting phase becomes quantum disordered. (See Figure 
(Id)). In this region, the properties of the model are better described by first diagonalizing 
the kinetic energy term of the Hamitonian. Both the triplet magnetic excitations and the 
doublet pairing excitations have finite energy gaps. The transition from the quantum dis- 
ordered phase into the antiferromagnetic phase is second order, the doublet pairing mode 
of the quantum disordered phase naturally continues into the 7r doublet mode in the anti- 
ferromagnetic phase, while the triplet magnetic excitation of the quantum disordered phase 
becomes the two gapless Goldstone modes of the antiferromagnetic phase. The transition 
from the quantum disordered phase into the superconducting phase is second order as well, 
the triplet magnetic excitation of the quantum disordered phase naturally continues into the 
7r triplet mode of the superconducting phase,0 while the doublet pairing mode becomes the 
gapless Goldstone mode of the superconductor. 



2 I would like to thank Prof. T.M. Rice for a stimulating discussion which lead me to this 
identification. 
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The SO (5) quantum disordered phase differ from the 50(3) quantum disordered phase 
of the quantum antiferromagnet [Q and the quantum XY disordered phase [I0| since it 



has both a massive triplet magnetic excitation and a massive doublet of pairing excitation. 
It is tempting to identify it with the resonating- valence-bond (RVB) phase introduced by 
Anderson. In fact, the type (Id) phase diagram looks very similar to the original RVB idea 
proposed by Anderson fl|]. The basic idea is that with low dimensionality and increased dop- 
ing, quantum fluctuations will first "melt" the antiferromagnetic state to form a resonating 
valence bond liquid of spin singlet pairs, and these singlet pairs further condense into a su- 
perfluid state. However, the 5*0(5) quantum disordered phase differ from Anderson's RVB 
phase in terms of the excitation spectrum. In the RVB phase, spinons and holons are the 
elementary excitations, while in the .SO (5) quantum disordered phase, even the electrons 
are "bound" into the collective coordinates. In some sense, the 50(5) quantum disordered 
phase describe a "incompressible liquid" of singlet pairs, similar to the Laughlin liquid in 
the fractional quantum Hall effect. 

We therefore see that the global features of the phase diagram deduced from the 50(5) 
quantum nonlinear o model agree reasonablly well with the general topology of the exper- 
imentally observed phase diagram of high T c superconductors. With this model, we can 
understand the depression of the Neel and the superconducting temperature in terms of a 
common origin, namely a region in the phase diagram where the 50(5) fluctuations are 
maximal. This theory makes precise predictions about the scaling behavior in the crossover 
region. These predictions can be tested experimentally in sufficiently clean systems where 
a direct first order transition becomes possible. I believe that phase diagrams of the type 
(la) or (lc) are reasonablly close to the reality in the high T c superconductors. Type (lb) 
may be useful to understand the heavy fermion materials, while type (Id) could be realized 
in the quasi one dimensional ladder systems where quantum fluctuations are strong. 

Theory of the spin gap, the superconducting transition and their relationship 
One of the most interesting phenomenon of the underdoped high T c materials is the so called 
spin gap. Above a certain temperature T s , the magnetic correlation length increases with 
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decreasing temperature, while it saturates below T s > T c . This phenomenon has only been 
observed in the YBCO family of the high T c superconductors. The spin gap phenomenon 
has been observed in two different kinds of nuclear spin relaxation processes, the 1/Ti and 
I/T2G measurements respectively. Through a series of insightful phenomeno logical analysis 
of the NMR experiments, Sokol and Pines |57 concluded that this two experiments can 



be linked through a simple z = 1 dynamical scaling hypothesis. Since our model is Lorentz 
invariant, their hypothesis is valid here. Therefore, we shall only address here the I/T2G 
experiment, which meausures the static spin correlation length. 

The phenomenon of spin gap has a natural explanation within the current theoretical 
model. There are three different temperature scales in the current model. When the temper- 
ature is lowered below T MF , the superspin vector acquires a finite magnitude. However, for 
a finite range of temperature below T MF , say T s < T < T MF , (see Fig. 1), the temperature 
scales is still high enough so that the superspin does not "notice" the anisotropy is its ori- 
entational degrees of freedom. In this range of temperature, the model is essentially 5*0(5) 
symmetric, and the antiferromagnetic correlation length increases together with the pairing 
correlation length as temperature is lowered. For T < T s , the thermal energy becomes low 
compared to the anisotropy energy in the superspin space, and the superspin vector lies 
preferably in the superconducting plane. Therefore, below T s , the antiferromagnetic corre- 
lation length saturates to a finite value. Finally, at T c , the superspin vector picks a particular 
direction within the superconducting plane, therefore breaking the U(l) gauge symmetry, 
and the system becomes superconducting. A very similar picture applies to the antiferro- 
magnetic side of the transition. There is a "pairing gap" temperature T p , above which the 
pairing correlation length increases together with the antiferromagnetic correlation length, 
while it saturates below T p . 

Both T s and T c can be calculated quantitatively within the 5*0(5) nonlinear a model. In 
principle, one can perform the calculation for all four types of phase diagrams shown in Fig. 
1. Due to the space limitation, we will only show the calculation for a representative case, 
the case (lc), where the transition from antiferromagnetism to superconductivity occurs via 
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a zero temperature quantum critical point. Since the anisotropy is irrelevant near the critical 
point, we will take the model to be isotropic. The idea of the large N mean field theory is 
to implement the constraint in the partition function (|TJ|) through a Lagrangian multiplier 
A and determine the value of A through a self-consistency condition that n 2 = 1 is enforced 
on average. Skipping the intermediate steps ||, the resulting self-consistency condition is 
given by 

r' 1 d 3 q 1 



o (2tt) 3 VcY + d 2 ^' 1 tanh(^d 2 q 2 + d 1 ^' 1 ) 

f 1 ' 1 d 3 q 1 . , 

+3 / . * . 30 

^0 (2tt) 3 ^J c 2 q 2 + c 2^-2 + ^2 tanh (^ c 2 q 2 + ^-2 + ^2) 

Here a is a lattice cutoff, £ is the superconducting correlation length, u; 2 = (2/x) 2 — (2/i c ) 2 
is the energy of the 7r triplet mode in the superconducting state, £,T 2 = £ -2 + {ojq/c) 2 is the 



antiferromagnetic correlation length and c = y p/x 1S t ne sound or spin wave velocity. 

The quantum bicritical point is determined by the condition that at ujq = 0, £ diverges 
at zero temperature. This determines 

5 

VpX = ^ a ' 2 (31) 

to be the condition for the quantum bi-critical point. When u>o is greater than zero, £ 
diverges at a finite temperature T c . The relationship between these two quantities can be 
easily determined from equation (|30|) , and is found to be 

^(^ 2 -3/(x))^ 16.44 (32) 

where x = ca~ l /ojq = £ s /a is the antiferromagnetic length at zero temperature, and f(x) = 
^xy/x 2 + 1 — hln(x + \/x 2 + 1). We see here that up to a logarithmic term, co>o and kT c 
are simply proportional to each other. This fact can be understood without any explicit 
calculations. At u = 0, the correlation length at a finite temperature is simply given by 
the Lorentz transformation of the finite width f3 in the imaginary time direction, or the 
thermal de Broglie wave length hc/kT. At T = and finite ujq, the correlation length can be 
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estimated from power counting the scaling dimension of uo, which gives c/ujq. Equating these 
two length scales and neglecting the logarithmic corrections gives the estimate of kT c oc u . 

Let us now turn to the calculation of the spin gap temperature T s . From the form of the 
antiferromagnetic correlation length £~ 2 = £~ 2 + (ujq/c) 2 , (a very similar form was first used 



see that at high temperature, the second term is smaller than the first, therefore, £ s depends 
on temperature in the same way £ does, namely proportional to the thermal de Broglie wave 
length hc/kT. However, at low temperature, the second term dominates, and £ s saturates to 
a finite, temperature independent value. The spin gap temperature can therefore be defined 
by c 2 £~ 2 (T s ) = u 2 . Inserting this definition into equation (0), we obtain a explicit formula 
for the spin gap temperature 



ducting transition temperature T c , the spin gap temperature T s and the resonant neutron 
scattering peak ujq to each other. £ s is the zero temperature antiferromagnetic correlation 
length averaged over three space dimensions. Since there are some uncertainies about the 
correlation length in the c axis and the value of the microscopic cutoff of the a model, we 
shall take £ s /a as a fitting parameter. 

The relation between the energy of the resonant neutron scattering peak and T c has 
been measured recently in the Y Ba^Cu^O^x superconductors by Fong et al fl28f . It was 
found that Uq is linearly proportional to T c with a slope of about 5.6. If we take £ s /a to 
be 0.8 in equation ( |32"D and roughly doping independent in the doping regimes explored in 
experiment, our theory gives the same slope. Notice that since £ s /a is a three dimensional 
average, the actual antiferromagnetic correlation length in the two dimensional sheet could 
be much greater. 

Combining the result of this section and the previous section on the collective modes, 
we see that the current work offers a unique explanation of the recent neutron scattering 



by Millis, Monien and Pines in their phenomenological analysis of the NMR data), we 




(33) 



Equations (^) and ( p3|) relates three experimentally observable quantities, the supercon 
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experiment of the high T c superconductors in terms of the pseudo-Goldstone bosons asso- 
ciates with the spontaneous 5*0(5) symmetry breaking. Since these magnetic excitations 
occur only below the superconducting transition temperature, the interpretation in terms 
of the Goldstone bosons is most natural, and we see that the linear proportionality between 
T c and u is also naturally explained within the current theoretical model. A class of the- 
oretical models |27| , |39| -p],[7i^2l attempt to associate the neutron resonance peak with the 
opening up of a d wave pairing gap, and possible excitonic states inside this gap. At the 
time when this class of theories were proposed, only a AlmeV peak was observed in the 
optimally doped YBCO superconductor p5|-p7[|. However, recent photoemission ||3| and 



resonant neutron scattering @,[29[ experiments in the underdoped regime can not be recon- 
ciled with this class of theories. In the underdoped regime, the electronic energy gap opens 
up at a temperature much higher than T c . If the resonant neutron scattering peak arises 
simply from the removal of spectral weight due to the electronic gap opening, such a peak 
would be observed far above T c , in contradiction to the experimental observation f28| , |2"9"|| . In 
this class of theories, the energy of the neutron scattering peak would also be proportional 
to the electronic gap energy. As the doping level is decreased from the optimal value, the 
electronic energy gap actually increases slightly P3|. Therefore, this class of theories would 



predict that the energy of the neutron resonance increases with decreasing doping, also in 



contradiction with the recent experiment by Fong et al. [[28 . 

This value of £ s /a « 0.8 can now be used to estimate the spin gap temperature without 
any additional fitting parameter. From equation (|3^) and (j33j), we find T s /T c m 1.37. 



There are some experimental uncertainties about the value of T s as well. The most acurate 
measurement of T s is in the T c = 62K material Y Ba2Cu^OQ,^- The \jTiG measurement 
shows that the antiferromagnetic correlation length saturate at T s « lOOf^ ||45|| . This gives 
a T s /T c ratio of 1.61, about 15% higher than the theoretically predicted ratio. 

We would like to remark that these quantitative calculations are strictly speaking only 
valid near the critical doping, while most of the experiments quoted above are far from it. 
It would be highly desirable to have sufficiently clean systems close to the critical doping 
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where precise comparison with theory can be made. 

Attempt at a synthesis The most important message from this work is that anti- 
ferromagnetism and superconductivity are the two different sides of the same coin. They 
are intimately related by a SO (5) symmetry which determines the competition and the 
low energy dynamics of the high T c superconductors. The current theory offers a unified 
explanation of some of the most challenging problems of high T c superconductivity. 

This work answers a number of experimental puzzles. The mechanism of superconduc- 
tivity is devided into a "high energy" piece, identified with the pair binding, or the formation 
of the amplitude of the superspin and a "low energy" piece, involving a "superspin flop" 
mechanism which selects the orientation of the superspin, and resolves the competition be- 
tween superconductivity and antiferromagnetism. These two different energy scales have 
different doping dependence. With increasing doping, T MF decreases because the effective 
J decreases, but T c increases because both the thermal and the quantum SO (5) fluctua- 



tions decreases. Recent photoemission experiment indeed found that the electronic gap 
decreases with increasing T c , a fact inconsistent with the weak coupling BCS theory, but 
consistent with the current model. This work also answers the question of why T c is low 
compared with the energy scales of the pair formation. Emery and Kivelson JTTJ rightfully 
argues that this happens because the superfluid density is low, but did not offer a explana- 
tion of why this is so. The current work explains this fact in terms of increased quantum and 
classical SO (5) fluctuations near the isotropic point /i c , so that both the Neel temperature 
and the superconducting transition temperature are depressed in this transition region. The 
puzzeling spin gap phenomenon is naturally interpreted in terms of the competition between 
entropy and anisotropy energy in the superspin space. So far the most direct experimental 
evidence for the approximate 50(5) symmetry is the resonant neutron scattering peaks in 
the underdoped and optimally doped region. These modes are naturally interpreted in terms 
of the pseudo Goldstone bosons associated with the spontanuous 50(5) symmetry breaking. 
Current model predicts that their intensity vanish exactly at T c , the mode energy is tem- 
perature independent, and decreases with decreasing doping. All these facts are consistent 
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with the experimental findings. Notablly absent in our discussion is the transport properties 
of the high T c oxides. Below T MF , transport properties could be addressed by considering 
fermions coupled to the 5*0(5) nonlinear a model degrees of freedom. However, quantitative 
details have not yet been fully worked out. 

The current theory also unifies a number of seeeming divergent theoretical approaches to 
the high T c problem. It outlines a strategy to systematically extract the low energy content 
of the t — J model by constructing low energy field theory constrained by the microscopic 
symmetry. The 5*0(5) symmetry can be used as a basic principle to organize the various 
theoretical proposals. Most directly, it unifies the 50(3) nonlinear a model theory HHH 
with the U(l) nonlinear a model theory fl0| , |TT|j . The "spin fluctuation exchange" fl2,|8|| 



approach should be interpeted as a theory of T MF . the phase separation in the t — J model 
Q occurs because the "superspin flop" transition is first order, therefore giving rise to a 
"forbidden density region" . The "spin bag phase" H can emerge in the phase diagram as 
a result of increased n fluctuation, or more precisely when Xn > Xc- The 50(5) quantum 
disordered phase could in some sense be associated with the RVB phase in Anderson's ]]J 
original proposal, and occurs as a result of increased 50(5) quantum fluctuation. 

The 50(5) theory makes a number of new experimental predictions. The most direct 
prediction is that when the materials are sufficiently clean, there could be a direct first order 
transition between the antiferromagnetic and the superconducting states. Measurements 
close to the bicritical point can be compared quantitatively with the theoretical predictions 
and can serve as important test of the theory. There are other predictions, due to the 
space limitations I shall only outline the basic ideas, and publish the details elsewhere. 
In the conventional Landau-Ginzburg theory of superconductivity, the order parameter is 
constrained to lie in a plane. Near the core of a superconducting vortex, a mathematical 
singularity is unavoidable. However, in the superspin model, since the order parameter is five 
dimensional, the superspin could lie in the superconducting plane far away from a vortex, 
but flip into the antiferromagnetic sphere inside a vortex. Such a topological configuration 
is called a "meron" in the field theory literature, meaning half of a Skyrmion. Therefore, 
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the current theory would predict that the core of a vortex in underdoped superconductors 
are not filled with normal electrons but are antiferromagnetic instead. Such a effect could 
be observed in the dissipation dynamics, or more directly, by studying the form factors 
of the vortex lattice in neutron scattering experiments. The current theory also predicts 
that at the interface between a antiferromagnet and a superconductor, the massive n modes 
become gapless, in very much the same way the collective modes in a Laughlin liquid in the 
quantum Hall effect become gapless at the edge. Such effects could be observed in tunneling 
experiments. 

High T c superconductivity is one of the most complex and perplexing phenomena ever 
encountered in condensed matter physics. Obviously, if there is any hope in completely 
solving this problem, the solution has to be derived from a simple and unifying principle. The 
present work is a attempt towards a synthesis of the complex high T c phenomenology from a 
symmetry principle which unifies antiferromagnetism with superconductivity. Throughout 
the history of physics, symmetry principle has been a faithful companion in our quest of the 
deep unity and harmony of Nature. May this time honored, deeply aesthetic principle help 
us one more time, to unlock the secret of high T c superconductivity. 
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FIGURES 

FIG. 1. Global phase diagram of high T c superconductors in the absence of disorder. Tmf is the 
temperature below which electrons bind to form singlet pairs, in the current model, it corresponds 
a finite magnitude of the superspin. T/v and T c are the Neel and superconducting transition 
temperatures respectively. There are four possible types of transitions from the antiferromagnetic 
to the superconducting state, a) There is a direct first order transition which terminates at a 
bi-critical point T& c . The first order transition can be described as a "superspin flop" transition, b) 
There are two second order phase transitions with a intermediate "spin bag" phase. The four second 
order lines merge at a tetracritical point Tt c - c) There is a single second order phase transition 
at a quantum critical point, d) There are two second order quantum phase transitions with a 
intermediate quantum disordered phase. In each type of the phase diagrams, there are two crossover 
temperatures T s and T p . The "spin gap" temperature T s corresponds to the temperature below 
which the superspin lies within the the superconducting plane, and "pairing gap" temperature T p 
corresponds to the temperature below which the superspin lies within antiferromagnetic sphere. 
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